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I.  STATBMEMT  OF  TBE  PROBLEM 

There  have  been  three  prieary  objectives  o£  our  recent  contract  to  ARO: 

The  developeent  of  analytical  derivative  methods  for  coupled-cluster  (CC) 
and  nany-body  perturbation  theory  (MBPT)  vave functions.  Such  methods  are 
essential  in  the  location  of  minima  and  transition  states  on  potential  energy 
surfaces,  and  in  the  first-principle  prediction  of  vibrational  spectra. 

Applications  of  the  CC/MBPT  methods  and  associated  analytical  derivative 
methods  to  transient  molecules,  such  as  those  present  in  the  combustion  of 
rocket  fuels,  in  interstellar  space,  or  occur  as  reaction  intermediates. 

The  development  of  CC/MBPT  methods  for  large  molecule  applications. 

Ve  have  made  a  number  of  achievements  in  these  three  areas.  The  folloving 
highlights  several  of  then. 

II.  SUMMARY  OP  SOME  IMPORTANT  RESULTS 

A.  During  the  course  of  this  contract  ve  presented  the  first  report  of 
MBPT<3)  and  D-MBPT<4)  (i.e.  fourth-order  MBPT  limited  to  double 
excitations)  analytical  gradients  (Bl).*  These  vere  illustrated  in 
calculations  on  H2O  and  CH2.  In  our  earlier  vork  for  ARO  (1-4]'*^  the  D- 
NBPT(4)  method  subject  to  finite  nuclear  displacement  vas  found  to  be 
quite  accurate  in  structure  determinations  for  flame  species  like  HCO, 
H2CO,  B^CO,  CH2N,  among  others.  Hence,  the  analytical  equivalent  should 
be  an  inexpensive  but  quite  useful  level  of  approximation  (see  E  below). 

*B  references  refer  to  the  bibliography  of  papers  supported  by  ARO  on  page  20. 
'''Other  references  are  listed  on  page  24. 
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B.  Ve  also  reported  the  first  CCD  (coupled-cluster  double)  analytical 
gradients  (B3]t  and  illustrated  this  aethod  to  H2CO.  CC/HBPT  aethods 
scale  correctly  vith  aolecular  size  but  are  non-variational. 

Consequently,  the  siaplifications  that  occur  for  variational  aethods 
cannot  be  used,  requiring  that  soae  rather  different  aethods  be  develc  ed 
than  in  SCF,  MCSCF  and  Cl. 

C.  The  above  application  of  CC/MBPT  theory  eaployed  a  procedure  that  would 
require  roughly  the  equivalent  of  a  CC  or  HBPT  calculation  for  each 
distinct  vibrational-rotational  degree  of  freedoa.  This  is  not  a 
practical  way  to  solve  the  analytical  CC/MBPT  problea,  so  ve  presented  a 
very  general  foraulation  that  deaonstrates  that  all  energy  derivatives 
vith  respect  to  3N  nuclear  displaceaents  nay  be  conputed  vith  only  about 
twice  the  effort  required  by  the  CC/HBPT  energy  calculation  [B4].  Only 
this  slaplification  nakes  It  possible  to  use  the  very  accurate  CC/MBPT 
aethods  for  searching  energy  surfaces  for  nolecules  with  several  degrees 
of  freedoa. 

D.  Building  upon  this  developaent,  ve  presented  a  paper  vith  detailed 
equations  for  any  level  of  CC  theory  and  extracted  the  full  fourth-order 
(i.e.  SDTQ-MBPT(4))  gradient  fomula  as  a  special  case  (B8]. 

B.  This  foraulation  is  very  different  than  our  original  analytical  gradient 
aethod,  and  it  requires  substantial  conputational  effort  to  implement. 
Just  recently,  ve  have  achieved  its  implementation  at  the  MBPT(3)  and 
D-MBPT(4)  levels  (B14].  Since  both  are  recognized  as  a  very  high  level 
of  correlated  theory  for  aolecular  structures  and  transition  states, 
computer  codes  ve  have  written  should  be  extensively  applicable  to  a  vide 
range  of  chemically  interesting  problems. 


F.  In  another  totally  new  development,  ve  developed  the  first  methods  for 
the  analytical  evaluation  of  HBPT(2)  second  derivatives  IB5,B6].  No 
routinely  applicable  correlated  second  derivative  method  vas  available 
prior  to  this  work.  Initial  illustrations  were  to  H2CO  [B5]  and  N2H2 
IB6J. 

6.  In  a  serious  application  of  our  analytical  MBPT(2)  second  derivative 
methods  to  a  large  basis  set  study  of  the  interstellar  species  SiC2,  ve 
shoved  that  the  triangular  fora  of  this  molecule,  observed  experimentally 
by  Smalley,  vas  indeed  a  minimum  on  the  HBPT(2)  energy  surface  [B7]. 

This  conclusion  differs  from  that  of  Grev  and  Schaefer  [5]  vho  found  the 
triangular  form  to  correspond  to  a  transition  state  at  the  SCP  or  GV6 
level,  but  they  assumed  that  a  higher  level  of  theory  vould  find  the 
triangular  structure  as  a  minimum,  once  substantial  amounts  of 
correlation  were  included.  Via  MBPT(2),  ve  find  the  triangular  form  to 
be  a  true  mininua.  Ve  also  predict  the  experimental  vibrational 
frequencies.  These  are  found  to  be  in  excellent  agreement  vith  an 
experiment  (6]  vhich  ve  found  after  our  results  had  appeared  (See  Table 
1). 


H.  Since  analytical  energy  second  derivatives  vith  respect  to  nuclear 

2 

displacement  (i.e.  8  E/9X  9T.)  provides  the  harmonic  force  constants  for 

a  p 

molecules,  a  generalization  to  the  mixed  second  derivatives  relative  to  a 

2 

nuclear  displacement  and  an  electric  field,  (i.e.  3  E/3X  3F.)  provide  the 

OK  Z 

dipole  derivatives  required  to  obtain  the  intensity  of  a  vibrational 
transition.  We  have  recently  added  these  terms  to  provide  MBPT(2)  level 


-3- 


Table  1 

SiCj  DZ+P  RING  STRUCTURE 


SCP^®^ 

GVB^*^ 

MBPT(2)^**^ 

Expt . ® 

Expt.' 

(C-C) 

1.256 

1.267 

1.294 

1.268 

1.250 

R(Si-C) 

1.835 

1.821 

1.835 

1.837 

1.812 

C-Si-C 

40.0 

40.7 

41.3 

40.4 

40.4 

C-C  str. 

1976 

1892 

1720 

1742® 

Si-C  A^  str. 

851 

873 

839 

852* 

Ring  Open. 

328i 

1711 

183 

^.S.  Grev  and  H.F.  Schaefer,  J.  Chea.  Phys.  80,  3552  (1984). 

^G.  Fitzgerald,  S.J.  Cole  and  R.J.  Bartlett,  J.  Chea.  Phys.  85,  1701  (1986). 

^P.  Tahddeus,  S.E.  Cummins  and  R.A.  Linke,  Astrophysical  Jour.  283,  245  (1984). 
^D.L.  Hichalopoulos,  N.E.  Geusic,  P.R.R.  Langridge-Saith,  and  R.E.  Smalley,  J. 
Chea.  Phys.  80,  3556  (1984). 

^.A.  Shepard  and  V.R.H.  Graham,  J.  Chem.  Phys.  82,  4789  (1985). 
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IR  intensities  as  veil  [B12].  These  techniques  make  it  possible  to 
routinely  predict  IR  spectra  for  molecules  while  Including  much  of  the 
essential  electron  correlation  effects. 

I.  The  unknown,  transient  molecule  occupies  an  essential  place  in  the 

understanding  of  the  chemical  bond  and  as  a  reaction  intermediate,  since 
it  is  isoelectronic  with  ozone  (0^),  cyclopropane  and  propene  It 

may  exist  in  three  forms,  a  cyclic  ring,  triaziridine;  an  open  form, 
triimide,  and  a  triazene  structure  (see  Fig.  1).  None  have  ever  been 
observed.  Using  our  MBPT(2)  first  and  second  derivatives  we  have 
determined  that  all  three  forms  of  N^H^  are  local  minima  on  the  potential 
energy  surface  and  defined  this  detailed  structure.  Using  high-level 
MBPT(4)  and  CCSD-4-T(CCSD)  methods  we  predict  that  the  triazene  isomer  is 
the  lowest  in  energy,  being  more  stable  than  the  triimide  form  by  17.30 
and  triaziridine  by  41.24  (see  Table  2).  Furthermore,  we  have  predicted 
the  IR  spectra  of  the  three  isomers  (see  Figs.  2-7)  at  the  SCF  and 
the  correlated  MBPT(2)  level.  Note  the  dramatic  differences  between  the 
SCF  and  the  correlated  HBPT(2)  predictions,  particularly  for  triaziridine 
and  triimide.  The  differences  pertain  to  the  intensities  auid  the 
frequencies.  Since  we  know  the  error  in  the  frequencies  at  the  MBPT(2) 
level  is  typically  less  than  5)C  compared  to  10-20Z  for  SCF  predictions 
(see  Table  3),  we  believe  the  HBPT(2)  results  to  be  sufficiently  accurate 
to  be  used  as  a  guide  to  experimental  spectroscopists  to  identify  in 
interstellar  space  or  in  matrix  isolation. 

J.  In  two  other  studies  [B9,B11]  we  have  used  our  MBPT(2)  analytical  second 
derivative  methods  and  high  level  <'  calculations  to  further  characterize 
the  decomposition  of  B2Hg  ■*  2BH^.  A  particularly  interesting  conclusion 
from  this  study  is  that  the  accepted  experimental  IR  spectra  (7]  of  the 
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All  above  calculations  used  a  DZP  basis  set. 
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Table  3 


AVERAGE  PERCENTAGE  ERROR  FOR  VIBRATKWAL  FREQUENCIES 
RELATIVE  TO  EXPERIMENT  (ANHARMONIC) 

(DZP  BASIS  SETS) 


SCP  MBPT(2) 

HNO  20.5  4.7 
CFj  13.5  3.1 
NjHj  15.5  4.2 
SlCj  6.7  1.3 
CH3NO2  12.9  4.9 
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transient  molecule  is  erroneous  I  Since  BH^  can  only  be  made  by 
dissociation,  the  experiment  produced  BH^  from  the  decomposition  of 
BH^CO.  After  subtracting  all  peaks  associated  vith  CO,  B2Hg,  and  various 
impurities,  one  strong  peak  and  three  weak  IR  peaks  remained.  These  were 
ascribed  to  BH^.  Our  calculation  for  BH^  (Table  4)  gives  good  agreement 
in  terms  of  intensities  and  frequencies  vith  three  of  the  peaks,  but  the 
strong  peak  occurs  at  an  incorrect  frequency  compared  to  our  calculation. 
Furthermore,  our  prediction  of  its  intensity  shows  that  it  would  be  less 
than  one-tenth  as  intense  as  the  weak  peaks,  meaning  it  could  not  have 
been  observed  in  the  experiment.  This  study  attests  to  the  power  of 
modern-day  ab  initio  methods  to  offer  a  reliable  voice  in  resolving 
experimental  discrepancies.  In  Table  5  are  presented  MBPT(2)  frequencies 
and  intensities  for  ^2^6  ^2^6*  exceptional  agreement  we 

achieve. 

K.  Under  this  contract  we  have  also  devoted  substantial  effort  to  approaches 

for  providing  ^  initio  correlated  predictions  of  energies  and  properties 

for  large  molecules.  Since  many-body  methods  scale  properly  vith 

molecular  size,  unlike  Cl  methods,  they  offer  an  appropriate  vehicle  to 

study  correlation  effects  in  large  molecules.  The  basic  idea  of  such 

large  molecule  methods  is  to  require  some  form  of  electron  pair 

localization  to  make  most  of  the  molecular  integrals  involving  disjoint 

parts  of  a  molecule  vanish  |B2].  By  so  doing,  the  difficulty  of  an  ab 

initio  calculation  scales  roughly  as  the  square  of  the  number  of  bonds 
2  6 

(-M  ),  rather  than  as  -N  where  N»M  is  the  number  of  basis  functions. 

In  a  series  of  studies  on  prototype  systems  C2H^  and  C2H2  we  have 

shown  that  even  very  crude  selecti^'ns  of  localized  orbitals  could  be  used 
vith  CC  methods  to  obtain  the  same  accuracy  as  the  conventional  choice  of 
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TABLE  4 t  NBPT(2)  haraonle  vibrational  fraquanciaa,  and  infrarad 
abaorption  intanaitiaa,  calculatad  for  ^^BH^  using  tha  6-31G** 
basis  sat.  Also  includad  art  tha  obsarvad  fundanantals.  and 
intansitias,  Units  ara  tha  aaaa  as  in  tabla  4. 


Tabic  5.  ealealatsC  MCmi)  karMaie  (r«c«*aci*s.  m  ,  aaC  iafrara4  akaarcttaa  iataaaitiaa,  t  ,  Car 


a|  Takaa  (r#a  ra(.  |£|.  Iataaaitiaa  ata  «ivaa  aa:  a-aaCiaa,  a-attaaf,  aa-vary  attaaf. 

kl  taca  ia  iaakla-karaaalc  ayfcaaiaatiaa. 
at  Cat  akaaiaaC. 


SCF  orbitals  [82].  Furthermore,  various  extensions  of  these  techniques 
can  be  proposed  that  will  restrict  the  correlation  treatment  to  certain 
regions  of  space.  There  is  the  potential  for  an  important  development  in 
this  area  for  large  molecules  that  should  be  pursued  in  future  work. 

L.  An  important  problem  in  the  theory  of  point  mutations  is  the  tendency  for 
bases  of  DNA,  like  cytosine  to  tautomerize  (see  Fig.  8).  The  cyt*  form 
is  a  rare  tautomer  of  cytosine  which  will  not  replicate  properly,  often 
leading  to  tumor  formation.  Using  some  of  our  large  molecule  methods,  we 
have  studied  the  effects  of  electron  correlation  on  this  tautomerization 
using  models  of  cytosine,  adenine  and  uracil,  as  well  as  cytosine  itself 
[810].  Ve  find  that  uniformly,  the  effects  of  electron  correlation  are 
greater  for  the  amino  (cyt)  than  the  imino  (cyt*)  forms.  Similarly, 
correlation  favors  lactia  forms  of  uracil  over  lactam  forms.  Unlike 
uracil,  the  very  slight  energy  differences  between  the  two  forms  of 
cytosine  prohibit  any  definitive  conclusions  about  relative  stability  to 
be  drawn  without  extensive  calculations,  but  it  is  apparent  that  electron 
correlation  which  accounts  to  ~1000  kcal/mol  in  cytosine  will  be 
essential  in  understanding  the  frequency  of  tautomerization.  Additional 
studies  are  devoted  toward  predicting  IR  spectra  of  such  complex 
molecules. 
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Tautomeric  Pairs  of  Nuclei  Acid  Bases  of  DNA 


I 


Fig.  8 
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and  R.J.  Bartlett,  J.  Chem.  Phys.  85,  5143  (1986). 
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Bartlett  and  V.N.  Lipscomb,  Chem.  Phys;  Lett.,  in  press. 
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Stanton,  R.J.  Bartlett  and  V.N.  Lipscomb,  in:  Molecules  in  Physics, 
Chemistry  and  Biology,  Volume  3,  ed.  Jean  Maruani,  Reidel,  Dordrecht,  The 
Netherlands,  to  be  published  fall,  1988. 

12.  "A  Search  for  the  N-H-  Molecule,"  P-  Magers.  E.A.  Salter,  R.J.  Bartlett, 
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IV.  R.J.  BARTLETT  INVITED  RESEARCH  PRESENTATIONS  SUPPORTED  BY  ARO 

(1984-1987) 


April 

1987 

National  ACS  Meeting,  Denver,  CO.  "Application  of  New 
Methods  for  Correlated  Ab  Initio  Studies  of  Large 
Molecules  Using  an  FPS  Scientific  Computer." 

March 

1987 

Parr/Eliel  Symposium,  Chapel  Hill,  NC.  "The 

Correlation  Problem:  From  PPP  to  Coupled-Cluster 

Theory . " 

June 

1986 

Canadian  Symposium  on  Theoretical  Chemistry,  Toronto, 
Ontario,  Canada.  "Analytical  Derivative  Methods  in 
Coupled-Cluster  Theory. 

May 

1986 

Southeastern  Theoretical  Chemistry  Conference, 
Tallahassee,  FL.  "Comparison  of  Intermolecular  Force 
Theory  and  Supermolecule  Approximations  for  Nonbonded 
Interactions. " 

April 

1986 

Workshop  on  Interface  between  Electronic  Structure  and 
Dynamics,  Snowbird,  Utah.  "Analytical  Derivative 
Methods  and  Coupled  Cluster  Theory." 

June 

1985 

Symposium  on  Molecular  Spectroscopy,  Ohio  State 
University.  "Many-Body  Studies  of  Molecular  Force 
Constants,  Infra-red  Intensities,  and 

Polarizabilities. " 

April 

1985 

Fulbright  Distinguished  Lecturer,  University  of 
Arkansas,  Fayetteville,  AK.  Two  lectures. 

Feb. 

1985 

Mardi  Gras  Symposium,  New  Orleans,  LA.  "Radical 

Studies  with  Many-Body  Methods." 

Dec. 

1984 

Nobel  Laureate  Symposium  on  Applied  Quantum  Chemistry, 
Joint  Meeting,  Pacific  Basin  Chemical  Societies, 
Honolulu,  HI.  "Radical  Studies  with  Many-Body 

Methods." 

Oct. 

1984 

NSF  Workshop  on  Future  Directions  for  Supercomputer  Use 
in  Chemistry,  Evanston,  IL. 

West  Coast  Theoretical  Chemistry  Conference,  Los 

Alamos,  NM. 

April 

1984 

March 

1984 

International  Symposium  on  the  Quantum  Theory  of 

Matter,  Palm  Coast,  FL. 

The  following  co-authored  papers  supported  by  ARO  were  presented  by  members  of 
the  Bartlett  research  group. 


Fifth  American  Conference  on  Theoretical  Chemistry,  June  15-20,  1984 

William  Laidig  Analytical  Gradients  for  the  Coupled  Cluster  Method 


Sanibel  Symposia,  March  14-23,  1985 
Robert  Harrison 

George  Fitzgerald 


Coupled  Cluster  and  MBPT  Analytic  Energy 
Se' '’nd  Deriva^ive£■ 

Gradients  of  Coupled  Cluster  Energies 
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George  Fitzgerald  Analytic  Energy  Gradient  Methods  for 

General  Coupled  Cluster  Energies 

David  Magers  Isomers  and  Vibrational  Spectra  of 


Snowbird  Conference,  April  21-25,  1986 

Ludvik  Adamovicz  Coupled  Cluster  Studies  on  Excited  States 

of  Polar  Diatomic  Anions 


Southeastern  Theoretical  Chemistry  Association  (SETCA)  Conference 
- May  1^ - - - 

E.A.  Salter  Analytic  Evaluation  of  Infrared 

Intensities 

D.  Magers  Structure  and  Spectra  of  N~H>:  Does  It 

Exist?  ^ 


V.  PARTICIPATING  SCIENTIFIC  PERSONNEL 

Graduate  Students:  David  Magers,  E.  Alan  Salter 

(Ph.D.  degrees  expected  1987) 

David  Bernholdt 
(Ph.D.  expected  1991) 

Postdoctoral  Research  Associates:  George  Fitzgerald,  Robert  Harrison, 

Sam  Cole,  Ludvik  Adamovicz 
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